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The symmetrized Density-Matrix- Renormalization- Group (DMRG) method is used to study linear 
and nonlinear optical properties of Free base porphine and metallo-porphine. Long-range interacting 
model, namely, Pariser-Parr-Pople (PPP) model is employed to capture the quantum many body 
effect in these systems. The non-linear optical coefficients are computed within correction vector 
method. The computed singlet and triplet low-lying excited state energies and their charge densities 
are in excellent agreement with experimental as well as many other theoretical results. The rear- 
rangement of the charge density at carbon and nitrogen sites, on excitation, is discussed. From our 
bond order calculation, we conclude that porphine is well described by the 18-annulenic structure 
in the ground state and the molecule expands upon excitation. We have modelled the regular met- 
alloporphine by taking an effective electric field due to the metal ion and computed the excitation 
spectrum. Metalloporphines have D±h symmetry and hence have more degenerate excited states. 
The ground state of Metalloporphines show 20-annulenic structure, as the charge on the metal ion 
increases. The linear polarizability seems to increase with the charge initially and then saturates. 
The same trend is observed in third order polarizability coefficients. 
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I. INTRODUCTION 

The electronic structure of porphyrin and metallopor- 
phyrins has been the subject of extensive study both be- 
cause of their interesting optical properties and chemi- 
cal activity [l|. Metalloporphyrins are biologically im- 
portant molecules, particularly those containing Fe (in 
Haemoglobin) and Mg (in chlorophyll) ions. The car- 
bon atoms that form the conjugated back-bone undergo 
easy substitution leading to a large class of porphyrin 
system with interesting electronic and chemical proper- 
ties. The 7T— conjugation present in these molecules gives 
rise to large nonlinear optic (NLO) responses, which can 
be tuned by peripheral substitution. Porphyrins are po- 
tentially useful in applications such as 3-D optical mem- 
ory devices 0] and optical power limiting Q. Solid por- 
phyrins are porous and could find use as molecular sieves 
and shape selective catalyst |4|. Porphyrins are chemi- 
cally and photochemically stable and are being used in 
photo-dynamic therapy [5j and antiviral therapy. Oxo- 
Vanadium(IV) porphyrins are studied for their anti-HIV 
properties @. Efforts are being made to fabricate solar 
cells with porphyrin based materials Q , since porphyrin 
framework is found in chlorophylls. Closely related to 
porphyrins are the phthalocyanines (Pc), which have a 
central porphyrinic core. The Pc system has attracted a 
great deal of attention in recent years for both organic 
electronic and spintronic applications [jj-Hlj. They have 
the interesting property of being adsorbed on metal sur- 
faces in a flat orientation. Scanning tunneling micro- 
scopic studies have been able to measure spin and charge 
densities at various positions in the molecule. 

Most of the important properties of porphyrins are 



dictated by their electronic structure. Hence, electronic 
structure of porphyrins and metalloporphyrins have been 
studied by many groups over several decades. Early stud- 
ies by Goutermann employed a four orbital model to ex- 
plain the nature of the Q band (the low-frequency weak 
absorption band) and the Soret (B) band (the high fre- 
quency intense absorption band) observed in porphyrins 
12]. Weiss et al. used the mean field results of Pariser- 
Parr-Pople (PPP) model Hamiltonian to explain the na- 
ture of Q and B bands of porphyrins [l3J . Christoffersen 
et al. have carried out ah initio calculation to study the 
electronic spectra of porphyrins, metalloporphyrins and 
other substituted porphyrins [14] . Zerner et al. showed 
that the Soret band is degenerate with parallel polariza- 
tion with respect to Q band, using random phase ap- 
proximation within the intermediate neglect of differen- 
tial overlap (INDO) method [1^]. Using density func- 
tional theory Ruth et al. have studied the singlet and 
triplet spectra of Zn-porphyrin and related compounds 
and showed that the B band consists of multiple absorp- 
tion bands [161 ] . 

Several groups have computed the NLO properties of 
modified porphyrins. Priyadarshy et al. |17| obtained 
the first hyperpolarizability of porphyrin bridged donor- 
acceptor molecules, wherein, the donor and acceptor 
molecules are attached at the meso— position. Albert et 
al. [l8[ have shown that large 2 nd order NLO response 
coefficients can be obtained by attaching a donor group 
at ft position and an acceptor group at meso— position; 
this is expected to increase the difference in the dipole 
moment between ground and excited states, leading to 
large first order polarizability. Shirk et al. have stud- 
ied the optical limiting properties of lead phthalocyanine 



2 



19]. Two-photon absorption (TPA) cross section of por- 
phyrin and Zn-porphyrins have been studied by a three 
state model by Zhou et al. [13] ■ Effect of donor-acceptor 
strengths on TPA cross section of aggregates of asym- 
metric Zn-porphyrins have been studied by Ray et al. 

a 

Porphyrin molecule is a fairly large system and hence 
is not amenable to accurate ab-initio studies. However, 
most of the interesting properties of porphyrin are as- 
sociated with the 7T— system. Thus the electronic prop- 
erties of porphyrins can be modeled by employing the 
well known PPP model for the ir— electrons. The Por- 
phyrinic 7T— system consists of 24 orbitals in conjugation 
occupied by 26 electrons. Even though the Fock space of 
the system is finite (4 24 = 2 48 - 2.8 x 10 14 ), the large 
Hilbcrt space of porphine, prohibits doing an exact quan- 
tum many body calculation. Specializing to desired total 
spin and occupancy still leaves the Hilbert space dimen- 
sion to be very large (9.27 x 10 11 ). Approximate tech- 
niques such as restricted configuration interaction (CI) 
for a few low-lying states is not very reliable, since it re- 
lies on cancellation of errors, as in single CI. Thus it is 
important to bring to bear novel and accurate methods 
for important low-lying states as well as NLO response 
coefficient of the systems. 

In this paper, we have employed the Density Matrix 
Renormalization Group (DMRG) method to carry out 
reliable model many-body calculations of the ground and 
excited state properties. In the next section we give 
a brief introduction to the DMRG method. In section 
III, we discuss the model Hamiltonian and implementa- 
tion of the DMRG method to porphines. In section IV, 
we present our results on the properties of the low-lying 
states of porphine. In section V, we discuss the dynamic 
nonlinear optic coefficient computed for the system. 



II. INTRODUCTION TO THE DMRG METHOD 

The DMRG algorithm involves efficient and accurate 
ways of truncating the insignificant degrees of freedom 
from the Fock space of the system. The DMRG method 
proceeds as follows. Given a model many body Hamil- 
tonian H which we wish to solve for a system with total 
number of sites N. We start with a small lattice of 2L 
sites and obtain the Hamiltonian matrix of this lattice 
in the real space basis. We obtain by standard methods, 
the lowest or the desired eigenstates of the Hamiltonian 
matrix. From the eigenvector we construct the density 
matrix p, of the subsystems with L sites by integrating 
out the states of the remaining L sites. The matrix ele- 
ments of the density matrix, p s '^, of the L sites block in 
the k th eigenstate is given by, 

Ps',s' = ^2 C s,e C s>,e (!) 



where the k eigenstate of the 2L sites problem is ex- 
pressed as 

v^ = EE c ">>i e >' (2) 

s e 

with \s) and |e) being the Fock space basis states of the 
system and environment blocks respectively and C k e are 
the associated coefficients. The density matrix p k,L is 
diagonalized and eigenvectors corresponding to m domi- 
nant eigenvalues of p k ' L are used to span the Fock space 
of the system block. If the order of the matrix p kL is 
M x M, then the mx M matrix O formed from m eigen- 
vectors as columns of the O is used for transforming all 
the operators of the system block, L, from the full Fock 
space to the truncated density matrix eigenvector basis. 
Thus if A is an M x M matrix of an operator A defined 
over the M dimensional Fock space of the system block, 
then A = OAO^ gives the renormalized m x m matrix 
representation of A. By interchanging the environment 
and system blocks, we also obtain renormalized matrix 
representation of all the operators in the former environ- 
ment block. 

We now add two new sites to the 2L size system at a 
convenient part of the system, usually in the middle. Us- 
ing as basis functions, the direct product of the density 
matrix eigenvectors of the two blocks of the 2L system 
and the Fock states of the new sites, we can obtain a 
matrix representation of the 2L + 2 site Hamiltonian. 
From the desired eigenstate -0fe.2L+2 we obtain the den- 
sity matrix p k > L+1 for the two blocks of size (L + 1) of 
the augmented (2L+2) site system and proceed as before 
to obtain all the renormalized operators on the (L + 1) 
sites and iterate until the desired system size is reached. 
This algorithm of constructing the system is called the 
infinite DMRG algorithm. In this algorithm, from the 
eigenvectors at any system size, we can also obtain ex- 
pectation values of any operator of interest. If the op- 
erator consists of product of site operators of different 
sites in the same block, then we need to carry through 
matrix representation of the product by renormalizing it 
at each step and use the renormalized matrix to com- 
pute the desired expectation value. If on the other hand 
the operator is a product of site operators of different 
sites on different blocks, we can use the renormalized 
matrix representation of the individual site operators to 
compute the expectation values. It is possible to tar- 
get more than one low-lying state simultaneously in a 
DMRG procedure. This is carried out by computing the 
density matrix for each desired states separately and us- 
ing the dominant eigenvectors of an appropriately aver- 
aged density matrix of these states. The choice of com- 
mon DMRG basis for different states permits comput- 
ing matrix elements such as transition dipoles between 
states. To access higher excited states of interest such 
as the l^Bu states in polyenes, we need to employ sym- 
metrized DMRG procedure and compute the desired ex- 
cited state as a low-lying state in an appropriate sym- 
metry subspace. For example, we could access a triplet 
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state as the lowest energy excitation in a space of odd 
parity, other parity symmetry corresponds to the invari- 
ance of the Hamiltonian when the spin space is rotated 
by 7r around the Y-axis. 

For molecular systems the infinite DMRG technique 
is not sufficiently accurate since the density matrix of 
the system block does not correspond to the density ma- 
trix of the target superblock. This can can be remedied 
by employing the finite DMRG algorithm. In this tech- 
nique the infinite algorithm is employed till the desired 
superblock is attained. After this point, the system block 
is augmented by one site and the Hamiltonian matrix of 
the superblock is set up using the larger system block, 
two new site and an environment block with one fewer 
site generated in the infinite algorithm. From the desired 
eigenstate of the Hamiltonian matrix, the system block 
is again augmented by one site and the whole process 
is repeated until the environment block reduces to a sin- 
gle site. From here on, the environment block and system 
block are interchanged and the process continued until we 
reach the size of environment and system block achieved 
at the end of the infinite DMRG algorithm. This process 
of progressively increasing the size of the system block 
and decreasing the size of the environment block, back 
and forth is called finite DMRG algorithm (22j . 

The DMRG technique has been shown to be highly 
accurate for quasi-one dimensional systems and has been 
widely employed in the study of the strongly correlated 
model Hamiltonians such as the Hubbard models [23j j . 
the PPP models [24T - [27| and the Heisenberg spin models 
[28l - l3"H . The PPP model has long range interactions and 
topologically the model can be viewed as a higher dimen- 
sional model. However, the interaction terms in the PPP 
model are diagonal in the real space and are hence well 
accounted for in the DMRG method. 

In this work, we report our studies on free base 
porphine (FBP) and metallo-porphines obtained from 
DMRG studies of their electronic structure. Porphine 
is modeled by PPP model with Ohno parametrization 
for the long-range interactions. We have studied the op- 
tical excitation spectrum by using symmetrized DMRG 
method [HI]. We characterize optically excited states by 
studying various properties of the states, such as, charge 
densities and bond orders. We have also studied both 
linear and nonlinear optical properties of the system. 



III. METHODOLOGY 

Porphine contains 4 pyrrolic rings connected by 4 CH 
groups as shown in Fig. [1] Nitrogen in the pentagons are 
sp 2 hybridized. There are two types of nitrogen atoms 
in the conjugation, one is the aza nitrogen and the other 
is the pyrrole like nitrogen. Lone-pair of electrons on N 
are in the plane of the molecule in aza-ring systems and 
do not participate in conjugation, whereas, they are per- 
pendicular to the plane of the ring and are involved in 
7r-conjugation, in pyrrole rings. The singly occupied p z 




FIG. 1: General structure of FBP. Sites are numbered in bold 
and transfer integrals (in eV) are given in blue colour. La- 
belling of the sites meso, a and f} are also shown in the figure. 
The system has D2 symmetry and all other bond transfers are 
given by symmetry. 



orbital of the aza nitrogens are involved in 7r-conjugation. 
Each carbon atom in the conjugation contributes one 
electron to 7r-system. Thus, in free base porphine, 26 
7T— electrons and 24 orbitals are involved in conjugation. 
It has been established that PPP model is well suited for 
describing the electronic properties of conjugated organic 
molecules. The PPP model Hamiltonian is given by 

Hi = X/*4<A CT+ ]L *« a i,<A> + h.c (3) 

i.a <i,j>,a 

H 2 = J~] Uihi >a n i(T ' 

i 

H 3 = 2J Vij{hi - Zi)(hj - Zj) 

i>j 

Hppp = Hi + H2 + H3 

Here Hi corresponds to the one-electron part of the 
Hamiltonian with the orbital energy, e, of p z orbital and 
Uj, the transfer integral between the orbitals i and j on 
the bonded pair of atoms in conjugation. We discuss the 
parametrization involving N-orbitals later. We have used 
experimental geometry for our calculations [33]. The 
transfer integrals, are calculated using a linear ex- 
trapolation, i.e. 



t is the transfer integral (— 2AeV) corresponding to 
C - C bond length of 1.397 A in benzene and (Sr)ij 
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is (fij — 1.397) A where is the length of the bond 
between atoms i and j. 

The term H2 in the PPP Hamiltonian incorporates 
on-site electron-electron repulsion energy and this corre- 
sponds to Hubbard interaction. The Hubbard parameter 
Ui is the energy cost for creating double occupancy of 
the orbital i. For carbon we have taken the standard 
value of 11.26 eV. The terms in H3 correspond to inter- 
site coulombic interaction between electrons in orbitals i 
and j, and we introduce local chemical potential Zi which 
is the occupancy of the orbital i that leaves the i th site 
electronically neutral. The value of z for C is always 1 
while that of nitrogen depends upon the occupancy of 
the orbital involved in conjugation, for the aza-nitrogen 
which contributes one electron to the 7r— conjugation, z 
is 1 , while for the pyrrole nitrogen which contributes two 
electrons to the it— conjugation, z is 2; Vij is calculated 
using Ohno parametrization [34|. 

We have parametrized e and U for nitrogen atom to 
fit the experimental excitation gaps and corresponding 
transition dipole moments for singlets in free base por- 
phine. The parameters which give the best fit are as 
follows: for aza nitrogens, €n 5 = e N 20 = —3.20 eV 
and f/jv 5 = Un 20 — 12.34 eV; for pyrrole nitrogens, 
e N 12 = e Ni3 = —14.0 eV and Un 12 — Un 13 = 15.00 eV. 
The value of tcN is taken to be —2.3 eV and —2AeV 
for C — N single and double bonds respectively. Other 
transfer integrals calculated from Eq. 0] are shown in 
Fig. [T] The nitrogen parameters reflect the more com- 
pact 2p orbitals of nitrogen compared to 2p orbitals of 
carbon. 

We have employed the DMRG methods to study the 
ground state (gs) and low-lying states of porphine. The 
topology of the transfer terms in porphines is not strictly 
one dimensional. The porphine system can be built up 
from a 4-site ring, by adding two sites at a time, in many 
different ways. The scheme shown in Fig. [2] provides 
a highly accurate scheme for building up the porphine 
framework, besides retaining the symmetries of porphine 
at every stage of the DMRG implementation. We have 
also carried out finite DMRG sweeps to improve the ac- 
curacy. The accuracy of the DMRG scheme is bench- 
marked against Huckcl model results for the system. In 
the DMRG studies, we have retained 210 density matrix 
eigenvectors (DMEVs) and carried out two finite DMRG 
sweeps. The order of the Hamiltonian matrix is about 
45, 000 after using two symmetries, namely, spin invari- 
ance and C2 symmetry. The matrix is sparse and sym- 
metric. We obtain a few low-lying states of the matrix 
using Davidson algorithm. Various properties such as 
charge and spin densities and b.o.s in different eigenstates 
can be computed since the matrices of the creation op- 
erators, the occupation number operators of all the sites 
and the b.o operators of all the bonds are renormalizcd 
at every step. 

In the molecular systems with strong electron correla- 
tions, it is difficult to identify the desired excited states 
such as two photon states or optically allowed excited 
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FIG. 2: A highly accurate scheme for building the porphine 
structure for DMRG calculations. At every step of the DMRG 
algorithm, we add two new sites shown by filled circle. Pos- 
itive integers correspond to the sites of the left block and 
negative integers to the sites of the right block. 



states, unless symmetry is exploited. The Hamiltonian 
is spin conserving, however, exploiting the conservation 
of both S 2 and S z (S corresponds to total spin ) within 
a DMRG scheme is difficult. Instead we have used spin 
invariance symmetry and S z conservation, to reduce the 
dimensionalities of the spaces. Spin invariance symme- 
try corresponds to invariance of the Hamiltonian when 
all the spins are rotated by tt around the y-axis, within 
the M s = (M s is eigenvalue of S z operator) subspace. 
The spin invariance operator divides the basis into two 
subspaces, one of even total spin and another of odd total 
spin. The spatial symmetry we have employed is the C2 
symmetry, even though the molecule has D2 symmetry. 
The C2 axis chosen by us passes through the center of the 
molecule and is perpendicular to molecular plane. The 
spin invariance symmetry commutes with the C2 symme- 
try and the resulting group has four one dimensional rep- 
resentations. Optical transitions are from states which 
are even under C2 symmetry to those that are odd, and 
retain total spin. The triplet (S=l) state is targeted as 
the lowest state in the M s — 1 sector. Since we are 
dealing with a nonzero M s subspace, the spin invariance 
symmetry does not exist in this case and we only employ 
the spatial symmetry. Besides, for computing properties 
such as transition dipoles and dynamic linear and non- 
linear polarizabilities, we need to express states in the 
different symmetry subspace by using the same DMRG 
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basis states. In order to retain proper balance of various 
states, we construct an average density matrix, averaged 
over density matrices of several low-lying states in each 
subspace and employ dominant eigenvectors of the aver- 
aged density matrix. The dynamic nonlinear optics re- 
sponse coefficient are computed using correction vector 
(CV) technique [35|. 

Conventional frequency dependent linear and nonlin- 
ear optical coefficient calculations rely on the knowledge 
of transition dipole moments between the gs and all ex- 
cited states. For exact calculation of NLO properties, 
the sum over states (SOS) method is most widely used. 
Generally, between 30 to 100 excited states, computed 
within a restricted CI scheme are used to obtain NLO 
coefficients. Soos and Ramasesha [36| introduced the CV 
method for calculating NLO coefficients within an exact 
diagonalization scheme which incorporates the contribu- 
tion of all higher excited states, without explicitly obtain- 
ing the excited state eigenf unctions. In this method, a 
correction vector (CV) is obtained using the gs wavefunc- 
tion and the Hamiltonian matrix. CV method involves 
computation of the first and second order correction vec- 
tors (j)^{ui\) and (jy?) (u)x,u)2), which are defined as, 



(H — Eq + fiwi + ir)|^ (1) ( Wl ) >= p,i\G > (5) 



(H-£ G + ^ 2 + zr)|<4 2) (wi.wa) >=Ajl4 1} M > ( 6 ) 

where H is the Hamiltonian matrix in the chosen many- 
body basis, Eq is the gs energy, u>i, 0J2 are the excitation 
frequencies and fii is the i th component of the dipole 
displacement operator, (/Zj =/Ij— < G\fli\G >), and T 
is the average lifetime of the excited states. It can be 
shown that (j)^' (ujx) and <j>^ (a>i, UI2) when expressed in 
the basis of the eigenstates of the Hamiltonian \R > are 
given by, 



("1) >= E 



<R\fLi\G> 



En — Er 



■ hbj\ + iT 



\R > (7) 



TABLE I: Comparison of excitation gaps and transition dipole 
moments between DMRG and exact calculations, in the non- 
interacting limit. Gaps are given in eV and transition dipoles 
are in a.u. 



DMRG 


Exact 


m = 180 


m = 210 


Gap \fi x \ \(iy\ 


Gap \fx x \ \n v \ 


Gap \/j, x \ \ny\ 


1.67 0.17 2.74 
1.82 1.07 0.09 
1.96 0.45 2.38 
2.05 2.90 0.24 


1.65 0.22 2.73 
1.78 1.47 0.32 
1.95 0.65 2.36 
2.03 2.72 0.51 


1.50 0.44 2.76 
1.56 2.37 0.52 
1.80 2.25 0.44 
1.87 0.42 2.35 



iff 



(-u)i - uj 2 , -wi) >) (9) 



where Pijki generate all permutations of i), (uj\,j), 

(u>2,k), (ui3,l) leading to 24 terms for 7 (with oj a = 
lui + L02 + 0J3). The linear polarizability components are 
given by 

ay(o;) = (<0f ) HI^|G> + <^ 1) (-^)l^|G^O) 

The 4>iS are exact within the Hilbert space chosen for the 
Hamiltonian. The linear algebraic equations that can be 
solved efficiently by a small matrix algorithm developed 
by Ramasesha [371 ]. 

Ramasesha et al. incorporated the CV technique in 
the DMRG method and have shown it to be robust [38|. 
In the CV-DMRG procedure, the average density matrix 
is constructed from a weighted average of the gs and the 
density matrix constructed from the correction vector. 
We will concentrate on the spatially averaged values of 
first and third order optical response. The expression for 
these are given by, 

3 1 

i=l ^ 
3 1 

= E T7^»3J + ^w) ( 12 ) 



EE 

S R 
IS > 



< S\p,j\R >< R\fii\G > 



(E R -E g + huji + iT)(E s — Eq + hoj 2 



ir) 

(8) 



Therefore |^>^(cji) > and |0y (wi,^) > can be readily 
used to compute linear and nonlinear frequency depen- 
dent polarizabilities. The third order NLO coefficients 
corresponding to sum frequency generation, in terms of 
these correction vectors, are given by, 



where {i}s are Cartesian indices (x,y,z). 

We have compared our DMRG results for the non- 
interacting Hamiltonian with exact Hiickel MO results 
for porphine. In the above calculations, at each itera- 
tion, the density matrix is constructed as the weighted 
average of the density matrix of five lowest eigenstates. 
Table HI compares the optical gaps and transition dipole 
moments for two different DMRG cut-offs, to, from finite 
DMRG calculations with exact results. We observe that 
to = 210 gives quite accurate results and keeping higher 
to increases the demand on computational resources with 
only marginal improvement in the results. DMRG results 
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slightly overestimate the higher excited states (« 0.15 - 
0.20 eV). Due to near degeneracy in the excitation lev- 
els, the components of transition dipole moments do not 
agree very well with Hiickel MO calculations although 
the magnitudes are in good agreement and the qualita- 
tive trend is maintained. All the studies that we report 
in this paper are based on retaining 210 density matrix 
eigenvectors and two sweeps of the finite algorithm. 



IV. RESULTS AND DISCUSSION 

This section is divided into three subsections. In the 
next subsection, we discuss results of our study on FBP 
followed by the metallo-porphine. In the last section, 
linear and nonlinear optic properties are discussed. 



A. Free base porphine 

We have obtained six lowest lying states in A and B 
subspaces for singlets as well as triplets of FBP using 
symmetrized DMRG method. In Table HTl we present the 
energies of singlets and in Table Hill that of triplets. The 
ground state energy is set to zero, so the energies re- 
ported are the excitation gaps from the gs. We have also 
presented experimental gaps and absorption intensities, 
wherever available. DMRG results correspond to isolated 
molecule calculations in the gs geometry. The experimen- 
tal optical gaps are usually red shifted by about 0.5 eV 
from the gas phase values. The first peak observed in the 
optical spectra is at 1.98 eV while that calculated is at 
1.66 eV. If we include the red shift of the theoretical gap, 
the agreement is off by about 0.8 eV. This suggests that 
the first peak in the spectra may be due to n — tt* tran- 
sition, and our model being a purely 7r-electron model 
cannot account for such a state. The 1.66 eV transition 
in the solid state would be shifted to about 1 eV cor- 
responding to approximately 1200 nm, and will not be 
observed in an optical spectra. The 2.86 eV absorption 
to the 2 1 B state is X-polarized (see Table HI)) while the 
nearly degenerate 3 X B and 4 1 B states at 3.89 eV and 
3.91 eV peaks have X and Y polarizations respectively. 
It should be noted that the low-temperature spectra of 
the free base porphine shows a split in the B band with 
equal intensity and separated by 0.03 eV (40| . Other the- 
oretical studies also show nearly degenerate levels at this 
energy. The 4.18 eV band {5 l B) corresponds to the iV 
band while the 4.83 eV absorption to 6 1 £? corresponds to 
the weak L band transitions observed in FBP. 

In Table Hm we give energies of several low-lying triplet 
states. The lowest singlet-triplet gap in the system is 1.42 
eV from our calculations, experimental value is 1.58 eV, 
which is in fair agreement. The lowest triplet-triplet ab- 
sorption, from our calculations is found at 1.67 eV, while 
experimentally two very close triplet-triplet absorption 
peaks are found at 1.58 eV and 1.65 eV [H, El). Our 
calculations show weak absorption peaks at 2.09 eV and 



TABLE II: Excitation gaps of six low-lying excited states and 
corresponding transition dipole moments compared with the 
experimental results for optically allowed states. Calculated 
oscillator strengths are given in parenthesis. Experimental 
results are from reference [41| . The experimental oscillator 
strengths are normalized with respect to most intense ab- 
sorption. Sum of calculated intensity (fj^+fiy) for the state 
3 and 4 is taken to be unity as they are nearly degenerate 
and all others are normalized with respect to it. Transition 
dipole between gs and all excited A states strictly vanishes by 
symmetry. 



State Optical gap 


Im*I 




Observed 


label 


(in eV) 


(in a.u 


) (in a.u.] 


value 


1 L B 


1.66 (0.16) 


0.05 


0.30 


1.98 (0.01) 


2 1 A 


2.94 








2 l B 


2.86 (0.26) 


0.39 


0.05 


2.42 (0.05) 


3 1 A 


3.40 








4 1 A 


3.53 








^B 


3.89 (1.00) 


0.48 


0.06 


3.33 (1.00) 


4}B 


3.91 (1.00) 


0.10 


0.59 




h x A 


4.01 








5 X B 


4.18 (0.21) 


0.04 


0.35 


3.65 (0.87) 


6 l A 


4.33 








6 1 B 


4.83 (0.01) 


0.00 


0.07 


4.25 (0.87) 



TABLE III: Excitation gaps of low-lying excited states and 
transition dipole moments in triplet manifolds. Note that all 
the excitations are from B to A space. Experimental results 
are from references [? 1 and 1431. 



State Triplet enerj 


;ies T-T gap 




My 


Observed 


label 


(in eV) 


(in eV) 


(in a.u. 


(in a.u.) 


value 


l 3 B 


1.42 










2 3 B 


2.45 


1.03 








3 3 B 


2.93 


1.51 








l 3 A 


3.09 


1.67 


0.00 


0.18 


1.58 


2 3 A 


3.51 


2.09 


0.03 


0.02 


1.65 


4 3 B 


3.55 


2.13 








3 3 A 


3.55 


2.13 


0.05 


0.08 




4 3 A 


4.26 


2.84 


0.00 


0.29 


2.82 / 2.96 


5 3 B 


4.28 


2.86 








5 3 A 


4.83 


3.41 


0.02 


0.12 


3.23 



2.13 eV and a much stronger absorption is predicted at 
2.84 eV. Experimentally two absorption peaks are seen at 
2.82 and 2.96 eV, close to the theoretical values. Theory 
also predicts a triplet absorption at 3.41 eV while exper- 
imentally a triplet absorption is seen at 3.23 eV. Overall, 
we find that the DMRG results for the 7r-electron model 
are quite consistent with experimental observations. 



1. Electron density and bond order analysis 

We have analyzed the eigenstates by computing charge 
densities and b.o.s in various states. In the gs, all /3 car- 
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FIG. 3: Electron density and b.o. for the gs of porphine. 
Charge densities are shown on the left part of the structure 
while b.o.s are shown on the right part. The site indices and 
bond indices are given in red and blue, respectively, for half 
of the system. The D2 symmetry gives value for remaining 
sites and bonds of the molecule. 

bon atoms have slightly more than one electron. The 
electron density at a carbon atoms in aza ring is 0.89 
and in pyrrole ring it is 1.08 and have opposite acceptor 
(donor) tendencies. At meso— carbon, the charge density 
is 0.98. This is consistent with the earlier semi-empirical 
calculations which show that the C meS o is electron defi- 
cient and Cp is slightly electron rich [18[. Nitrogen in 
the aza ring has 1.41 electrons and in pyrrole ring it has 
1.63 electrons. These results are in agreement with the 
earlier theoretical work 13]. The gs charge density and 
b.o. are given in Fig. [3J It is interesting to note that 
although aza ring nitrogen contributes only one electron 
to 7r -conjugation the electron density is 1.41. 

For excited states which are Y-polarized, electron den- 
sity is transferred from nitrogen in aza ring to the neigh- 
boring carbon atoms while for the the X-polarized excited 
state the electron is transferred from the carbon atom 
to the N atom. Except for the fourth and fifth excited 
states, for which there is considerable transition dipole 
moment, the change in electron density on C a and Cp 
sites are negligible. In the fourth and fifth excited states, 
there is considerable electron density redistribution be- 
tween carbon and nitrogen atoms. The electron density 
at meso sites remain almost the same in all the excited 
states. 

In pyrrole rings the opposite trend of electron trans- 
fer is observed except for the sixth excited state. The 



TABLE IV: Electron density p for gs and difference of electron 
density Sp for the optically allowed states w.r.t. gs. 



Site 


P 


Sp 


gs (A ) 


B\ B2 B3 B4 B$ Be 


Cp (aza) 
C/3 (aza) 
C a (aza) 
C a (aza) 
N (aza) 


1.01 
1.01 
0.89 
0.89 
1.41 


0.01 0.00 -0.00 -0.03 -0.04 -0.00 
0.00 -0.00 -0.01 -0.02 -0.05 -0.01 
-0.01 0.01 -0.01 0.04 0.03 0.02 
-0.01 0.02 0.01 0.05 0.05 0.01 
-0.02 0.00 0.02 -0.11 -0.08 -0.03 


Cmeso 
Cmeso 


0.98 
0.98 


0.01 0.00 0.03 0.03 0.03 0.01 
0.00 -0.00 0.01 0.02 0.03 0.02 


Ca (pyr) 
Ca (pyr) 

C/3 (pyr) 

C/3 (pyr) 
N (pyr) 


1.08 
1.08 
1.02 
1.02 
1.63 


-0.02 0.04 0.02 -0.03 -0.02 0.04 
-0.01 0.03 0.03 -0.03 0.01 0.02 
-0.00 0.04 -0.01 0.00 0.03 -0.00 
0.01 0.04 -0.05 0.03 0.02 0.03 
0.04 -0.19 -0.05 0.06 -0.00 -0.11 



TABLE V: Bond orders of gs and difference in b.o. for six 
low-lying states w.r.t gs, in the singlet manifold are given. 
Symmetry space of states are given in parentheses. 



Bond 


b.o 
G.S (1A) 


A{b.o) 
(1**)(1B) 


A{b.o) 
2 nd (2A) 


A(b.o) 

y d (2B) 


A{b.o) 
4 th (3A) 


A(b.o) 

5 th (4A) 


A(b.o) 
6 th (3B) 


cv 


- Cp (aza) 


0.81 


-0.01 


-0.01 


-0.03 


-0.15 


-0.16 


-0.04 


Cp- 


- C a (aza) 


0.39 


0.00 


0.01 


0.02 


0.10 


0.08 


0.03 


c,- 


- C a (aza) 


0.38 


0.00 


0.01 


0.01 


0.11 


0.08 


0.02 




- N (aza) 


0.54 


-0.01 


-0.02 


-0.03 


-0.06 


-0.07 


0.00 


c a 


- N (aza) 


0.52 


-0.01 


-0.01 


0.00 


-0.02 


-0.05 


-0.03 


c c 


- C m€ so 


0.59 


-0.03 


-0.01 


-0.02 


-0.10 


-0.03 


-0.01 


Cc 


- Cme.o 


0.58 


-0.03 


0.01 


0.02 


-0.06 


-0.02 


-0.04 


Cc 


~ Cmeso 


0.55 


-0.01 


-0.02 


-0.02 


0.02 


-0.04 


-0.03 


C c 


~ Cmeso 


0.56 


-0.01 


-0.03 


-0.06 


0.00 


-0.06 


-0.01 


C a - 


- Cp (pyr) 


0.57 


-0.01 


0.04 


0.04 


-0.02 


0.06 


0.02 


Cc - 


- Cp (pyr) 


0.55 


-0.01 


0.04 


0.05 


0.00 


0.05 


0.01 


<v 


- Cp (pyr) 


0.58 


-0.04 


-0.12 


-0.19 


-0.07 


-0.14 


-0.18 


C a 


-N (pyr) 


0.37 


-0.02 


-0.06 


-0.04 


-0.05 


-0.06 


-0.11 


Cc, 


-N (pyr) 


0.38 


-0.03 


-0.03 


-0.02 


-0.06 


-0.04 


-0.09 



electron is transferred to N when the polarization is in 
Y direction and is transferred from N, when it is in X 
direction. The magnitude of electron distribution is max- 
imum in second excited state. In the fifth excited state, 
there is hardly any change in the electron distribution at 
iV atom, but the charge is transferred between C a and 
Cp atoms. Deviation in electron density from the gs are 
given in Table IIVI 

It is known from the literature that the porphyrin can 
have two conjugation pathways: an inner 18-annulene 
pathway and an outer 20-annulene pathway (Fig. |4j. 
From our b.o. calculations we find that b.o.s are more 
uniform (0.55 ±0.03) along the 18 annulene pathway. We 
should expect porphine to behave more as an 18-annulene 
system than a 20-annulene system. The two Cp — Cp 
bonds have much higher bond order (0.81) and hence the 
20 annulene description of porphines can be discarded. 
This structure is in confirmity with X-ray structure 13311 
as well as other theoretical work by Weiss et al. |13j . 
From Fig. 21 we see that pyrrole N is not strongly in- 
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(a) (b) 

FIG. 4: Two possible equilibrium geometries in gs : (a) 18- 
sites annulenic structure and (b) 20-sites annulenic structure, 
shown by bold lines. 




FIG. 5: Electron density, spin density and b.o.s for the lowest 
triplet state of porphine. Spin and charge densities are shown 
on the left part of the structure while b.o.s are shown on the 
right part. The D2 symmetry gives value for other half of the 
molecule. 



volved in 7r-conjugation, Cp — Cp b.o. is small and hence 
we expect bond length of C^j — C^ in pyrrole rings to be 
large. In aza rings, N is involved in the ir— conjugation 
path and a reverse trend in b.o. is observed, namely, 
C/3 — C^ b.o., is large and we should expect a shorter 
Cfl — C/3 bond. The C Q — bond length is increased 
[44j in aza rings. 

The devitation from the g.s. bond order is given in 
Table |Vj In the first excited state, all the b.o.s become 



smaller compared to their values in gs geometry. The 
change in b.o. is only marginal in aza rings, while in 
pyrrole rings 0,3 — 0,3 bonds become weaker. All the 
three C — C b.o.s are uniform. The C — N bonds also 
become weaker. The meso— bonds connected to aza rings 
become weaker than those connected to pyrrole rings. 
Hence we expect the equilibrium geometry in the first 
excited state to have expansion of the whole 7r-system. 

In the second excited state, C/3 — C/3 bond of aza rings 
weakens and C a — C^ bonds become stronger. Meso 
bonds connected to aza rings become weaker and the rest 
become stronger compared to the 1 st excited state. The 
pyrrole ring becomes more pyrrolic, i.e. the C^ — C^ b.o. 
reduces considerably and the C Q — Cp bond becomes a 
strong bond while the C — N bonds get weaker. This 
bond pattern is in confirmity with the charge distribu- 
tion in the second excited state. The charge is depleted 
from the N and accumulated on C atoms. The b.o. pat- 
tern of third excited state is similar to that of second 
excited state. 

In the fourth excited state, Cp — Cp bonds of aza rings 
become very weak. The C a — Cp and the C — N bonds 
have almost the same b.o. (0.49 ± 0.01). This can again 
be attributed to the fact that in the fourth excited state, 
there is charge accumulation on C Q atoms and depletion 
at N atom. The bond alternation of the meso— bonds 
reduce. In pyrrole ring, all the three C — C bonds are al- 
most uniform. Note that in both the first excited and 4 th 
excited states, there is charge accumulation at pyrrole N 
and all the three C — C bonds have uniform bond order. 

In the fifth excited state, the aza b.o. and charge dis- 
tribution pattern are same as in fourth excited state. The 
charge density pattern at pyrrole ring is similar to that 
of the third excited state and hence the b.o. pattern 
is similar to that of third excited state, namely alter- 
nate single-double-single bond for three C — C bonds and 
weaker C — N bonds. 

In the sixth excited state, the charge depletion at aza 
N is not significant and the b.o. is slightly alternating 
for C a — Cp and C — N bonds. The C/3 — C^ bond is more 
like a double bond. The pyrrole ring is less pyrrole like 
and the C — N bonds become weak. 

In aza rings, when the charge is taken away from N, 
the b.o.s of C a — Gp bonds and C — N bonds would be 
uniform. On the other hand, when there is excess charge 
at N, the bonds alternate. In pyrrole rings, when there is 
charge accumulation at N, the three C—C bonds become 
almost uniform and when the charge is taken away from 
N, C a — C/3 bonds become stronger and C^ — Gp bonds 
become weaker. In all the excited states, either the aza 
C — N bonds have weak bond alternation or the pyrrole 
C - N bonds. 

Spin densities, charge densities and b.o.s are shown 
in Fig. [SJ in the lowest triplet state of FBP. The elec- 
tron density on Cp and C Q carbon atoms are compa- 
rable to that in gs while it decreases on the nitrogen 
atom of aza-rings and remains constant for N atom on 
pyrrole-rings, compared to gs. At C meS o sites, electron 
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density decreases in the lowest triplet state. The val- 
ues of Cp — Cp, C a — Cp and C a — C mcso b.o.s in the 
aza-ring are almost same as in gs while C Q -C mG so and 
C a — Cp b.o.s of pyrrole-ring increase while for Cp — Cp 
it decreases. C — TV b.o.s in aza-rings decrease but in- 
crease in pyrrole-rings. Spin density values on Cp atoms 
in aza-ring are small but are relatively higher (0.045) 
in pyrrole-rings. Spin density values are very high on 
Cmcso sites (0.16) and nitrogen atom (0.12) in the aza- 
rings while this value is very small on the nitrogen atom 
of pyrrole-rings. 



B. Metallo-porphines 

Metallic porphines are biologically important 
molecules. There are two categories of metallic 
porphines - regular and irregular porphines. In the 
case of regular porphines, the metals have closed shells 
and in the latter, they have partially hlled shells. Our 
study concentrates on regular metallic porphines. We 
have placed the metal ion at the center of the porphine 
molecule and it acts only as a source of electric field, i.e., 
we have not taken the overlap of metallic d— orbitals 
with porphine 7r— orbitals. The modified Hamiltonian 
to account for the effect of the metallic ion is given by, 
H = Hppp + qAi{hi - Zi)/r M , where q M is charge on 
the metal and (fij — z>) is effective charge operator at 
site i, is the distance between metal at the center 
and i th atom. In the presence of the metal ion, pyrrolic 
TV gets deprotonated and all the four five-membered 
rings become equivalent and hence the bond alternation 
observed in the neutral porphine molecule also vanishes. 
So, we have taken uniform transfer integrals for all 
the C — C bonds and a single value of ejv and Un for 
all nitrogens. The uniform ejv and Un used in the 
Hamiltonian are -7.8 eV and 12.34 eV respectively. We 
have varied the charge on the metal ion from 1 to 5 and 
studied the effect of oxidation state of the metal ion on 
porphine spectra, ground state charge distribution and 
nonlinear optical properties. 



1. Ground state charge density and bond orders 

Charge densities on porphine for different oxidation 
states of the central metal ion are given in Table I VII 
Charge density on the Cp atoms decreases while it in- 
creases on C a atoms. The charge density on TV atoms 
increases while those on C mcso atoms decrease. 

Bond orders in gs for metallo-porphine as a function 
of oxidation state on the metal ion are given in Table 
IVIII As in the FBP molecule, C^ — C^ bonds in the aza 
rings are the strongest bonds except for +5 oxidation 
state. As the oxidation state of the metal ion increases, 
Cp — Cp bond order decreases and the Cp — C a b.o.s in- 
crease. With the increase in qm, as the charge on N 
reaches its maximum value, the b.o. of all C — N bonds 



TABLE VI: Charge densities on unique sites of metallo- 
porphine for different oxidation states of the metal ion. 



Site 




Char 


je on 


metal site 







1 


2 


3 


4 5 


Cp (aza) 


1.01 


1.04 


0.99 


0.94 


0.87 0.80 


C a (aza) 


0.89 


0.79 


0.82 


0.88 


0.96 1.05 


N (aza) 


1.41 


1.73 


1.82 


1.88 


1.92 1.95 




0.98 


1.08 


1.04 


0.99 


0.93 0.88 


C a (pyr) 


1.08 


0.78 


0.82 


0.89 


0.97 1.06 


C/3 (pyr) 


1.02 


1.05 


0.99 


0.92 


0.85 0.77 


N (pyr) 


1.63 


1.79 


1.86 


1.91 


1.94 1.96 



TABLE VII: The variation in b.o.s of the porphine bonds 
with different oxidation states of the central metal ion. B.o.s 
of bonds for which the value is not quoted can be obtained 
by imposing D2 symmetry of the porphine. 



Bond 


Charge 


on metal site 


1 


2 


3 


4 


5 


Cp 


-Cp 


0.74 


0.72 


0.69 


0.62 


0.57 


c (i 


~ C a 


0.47 


0.50 


0.54 


0.59 


0.62 




Ca 


0.47 


0.51 


0.54 


0.59 


0.62 




-TV 


0.41 


0.34 


0.27 


0.22 


0.17 


Col 


-TV 


0.42 


0.34 


0.27 


0.22 


0.17 


Col 


Cmeso 


0.59 


0.60 


0.60 


0.59 


0.57 


Col 


Cmeso 


0.59 


0.60 


0.60 


0.59 


0.57 


Col 


Cmeso 


0.57 


0.58 


0.59 


0.59 


0.60 


Ca ~ 


Cmeso 


0.57 


0.58 


0.59 


0.59 


0.60 


C/3 


Ca 


0.55 


0.57 


0.58 


0.61 


0.62 


C/3 


C a 


0.54 


0.56 


0.58 


0.60 


0.61 


Col 


— TV 


0.35 


0.28 


0.22 


0.18 


0.14 


Col 


-TV 


0.33 


0.27 


0.22 


0.18 


0.14 


C/3 


-Cp 


0.61 


0.61 


0.59 


0.56 


0.52 



decrease sharply. Meso-bonds show minimum change; 
one set of b.o.s increases while the other set decreases, 
leading to small alternation in the meso— bonds. The 
asymmetry in the geometry of two rings reduces as the 
charge on the metal ion increases. The C — N bonds of 
aza rings have almost same b.o.s for all oxidation state 
while in pyrrole rings they are of same magnitude only 
for higher oxidation states. 
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Our calculations are in good agreement with earlier 
semi empirical results. Poveda et al. (45j found that 
Cp — Cp bonds are the strongest bonds and the — C a 
bonds are the weakest bonds in both Mg Octaethyl tetra 
phenyl porphyrin (OETPP) and ZnOETPP. For charge 
+2 we also find that the 0,3 — 0,3 bonds are the strongest 
bonds and C mcS o — C Q bonds are stronger than C Q — 
bonds. The C — N bonds become weaker in the presence 
of the metal ion. All the C - N b.o.s of FBP are al- 
ways higher than those in the metallo-porphine. Pyrrolic 
C — N bonds have a slight alternation in bond order pat- 
tern at lower oxidation states. FBP in the excited states 
show alternate b.o.s for C — N bonds either in the aza 
rings or in the pyrrole rings. While FBP has a 18- 
annulenic structure, metallo-porphine shows 20-annulene 
like structure when the charge on the metal ion increases. 



2. Optical properties 

The metallo-porphines have higher symmetry than 
FBP; former has D^h point group symmetry while the 
point group symmetry of the latter is D2h- Hence the 
spectra of metallo-porphines have more degenerate exci- 
tations than in FBP. 

The characteristic optical spectra of metallo-porphines 
consists of a Q-band lying between 500-600 nm. Q- 
bands have two nearly degenerate levels separated by 
1250 cm ■ It was originally identified as a vibronic pro- 
gression and assigned Qo->o and Qo->i transition [46j. 
B-band which is the most intense band lies between 380 
and 420 nm. Other L, M and N bands are comparatively 
weak and appear between 325 nm and 215 nm. 

The five lowest excitation gaps along with transition 
dipole moments are given in Table IVIIII for five differ- 
ent oxidation states of the metal ion. The first allowed 
state lies about « 1.70 eV above the gs. For <?m = +1 
and qM = +2 transition dipoles decrease with increas- 
ing charge. The lowest level has almost no intensity for 
dipole transition from gs for qM — +3 and the second 
level has significant intensity. We find that for metal 
ion oxidation states upto qM = +3, Q-bands are well 
separated. But with increase in qM, the separation be- 
tween Q-bands decreases. For qM > +3 on the metal ion, 
this separation is very small, and both states have large 
transition dipoles to the gs although with different polar- 
izations. The intensity for second Q-band transition for 
qM = +2 is weak and is consistent with experimental re- 
sults [47( | . Experiments show that the first excited state 
of metallo-porphines is slightly lower in energy than in 
FBP, irrespective of substitutions at the meso— positions 
as shown in reference (48j . There is only a small vari- 
ation (maximum of O.leV) in optical gap for different 
substituents. 

From our DMRG results, the strongly optically al- 
lowed band, B-band, lies between 3.8 ± 0.2eV above the 
gs for all values of qM- These levels show red shift of 
about 0.2eV as qM is increased from +1 to +5. For 



TABLE VIII: Excitation gaps for singlet-singlet transition of 
five lowest optically allowed states and transition dipole mo- 
ments for different metallic charges. 



qM 


Excited Optical gap 


M 


state 


(in eV) 


(in a.u.) 




IB 


1.70 


0.80 




2B 


2.04 


0.08 


1 


3B 


3.94 


0.66 




4B 


4.11 


2.30 




5B 


5.37 


0.26 




IB 


1.69 


0.46 




2B 


1.97 


0.16 


2 


3B 


4.00 


2.43 




4B 


4.22 


0.66 




5B 


5.36 


0.59 




IB 


1.70 


0.06 




2B 


1.84 


0.41 


3 


3B 


3.87 


2.53 




4B 


4.23 


0.58 




5B 


5.13 


0.55 




IB 


1.66 


0.66 




2B 


1.73 


0.71 


4 


3B 


3.59 


2.41 




4B 


4.36 


1.38 




5B 


5.20 


0.40 




IB 


1.65 


1.02 




2B 


1.68 


1.06 


5 


3B 


3.62 


2.26 




4B 


4.31 


1.32 




5B 


4.78 


0.10 



qM = +2, the experimental gaps are about 3.22eV for 
Zn and 3.18eV for Mg(etio)Porphine (42$. Time depen- 
dent density functional theory calculation of Minaev et 
al. has a blue shift of 0.23 [49|], while symmetry adapted 
cluster CI study of Hasegawa et al. has a blue shift of 
0.5 eV for the B band compared to the experimental re- 
sults [loj]. The blue shift in our DMRG studies is 0.6 eV. 
Complete active space second order perturbation theory 
(CASPT2) calculation of Rubio et al. underestimates the 
experimental excitation energies by 0.3 to 0.5el^ (5lj |. 

We observe a transition of medium intensity around 
4.2 ± 0.1 eV in all metallo-porphines, correspond- 
ing to experimental absorption peaks for iV-bands of 
Mg(etio)Porphine and TPPMgP observed at 3.85eV and 
S.Q7eV respectively. Our result overestimates this gap by 
fa O.SeV. There is also a weak L-band around 5.3±0.1eV. 
To the best of our knowledge, there are no experimental 
results for system with qM > 3, for regular metallopor- 
phines for comparison with DMRG results. 
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TABLE IX: Triplet state energies, T-T gaps and correspond- 
ing transition dipole moments of metallo-porphine for differ- 
ent oxidation state of metal ions. 



Qm 


state 


Triplet 


T-T gap 


M 




index 


energies (eV) 


(eV) 


(in a.u.) 




1 


3.31 


1.29 


0.30 


1 


2 


3.48 


1.46 


0.59 




3 


4.11 


2.09 


0.28 




4 


4.91 


2.89 


0.17 




1 


3.01 


1.27 


0.34 


2 


2 


3.33 


1.59 


0.29 




3 


3.84 


2.10 


0.25 




4 


4.41 


2.67 


0.19 




1 


2.99 


1.27 


0.19 


3 


2 


3.31 


1.59 


0.08 




3 


3.62 


1.90 


0.47 




4 


4.23 


2.51 


0.12 




1 


3.29 


1.69 


0.07 


4 


2 


3.45 


1.85 


0.09 




3 


3.99 


2.39 


0.48 




4 


4.11 


2.51 


0.24 




1 


2.76 


1.60 


0.69 


5 


2 


3.25 


2.09 


0.51 




3 


3.30 


2.14 


0.55 




4 


3.73 


2.57 


0.26 



3. Triplet- Triplet spectra 

In Table IIX1 triplet energies and triplet-triplet (T-T) 
gaps and corresponding transition dipole moments are 
given. Goutermann et al. [HiJ observed a weak T-T ab- 
sorption at lAleV and a high energy T-T absorption at 
1.58eT^ with high intensity in case of TPPZn. Sapunov 
et al. observed a peak at 1.53eV for the same system 
in dimethyl phthalate solvent [42j while Pileni et al. ob- 
served peaks at 1.48eV and 1.68eV in DODAC solution 
of TPPZn. Goutermann et al. also measured the phos- 
phorescence of (TPP)SnCl2, which exhibits similar be- 
havior as TPPZn. DMRG results for q M = 2 have the 
first T-T transition at 1.27eV and the second at 1.59eV. 
The first transition is red-shifted compared to the exper- 
imental value. In our calculation we observe that the 
intensity of the first peak is comparable to the second 
peak. Other higher excited states have relatively weak 
intensities. 

As the oxidation state of the metal ion increases, the 
first T-T gap remains almost a constant at 1.27eV up to 
qm = 3 and then it increases for charge qM = 4 and 5. 
All the T-T gaps for charge qm = +2 and qM = +3 are 
similar. For charge qM — 4, the first observable transition 
is at 2.39eV. For qM = 5, all the T-T transitions have 
significant transition dipoles. 



V. NONLINEAR OPTICAL PROPERTIES OF 
PORPHINES 

Hitherto, the dynamical NLO response of porphines 
has been studied using a sum-over-state (SOS) method 
in conjunction with restricted CI studies. As has been 
stated before, DMRG method is far superior to other 
many-body technique for quasi-one-dimensional systems. 
Besides, the DMRG method has been combined with the 
CV method to bypass computing only a part of the ex- 
cited state spectrum as required in the SOS method. The 
CV method incorporates the full excitation spectrum in 
the truncated DMRG basis, which in our scheme cor- 
responds to about 45,000 excited states for the DMRG 
cut-off of m = 210. 

We have calculated the linear and third order nonlinear 
optic coefficients corresponding to third harmonic gen- 
eration (THG). These r esponse coefficients have been 
calculated at three different frequencies. These compu- 
tational results have been obtained for both FBP and 
regular metalloporphines. 
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FIG. 6: Tumbling averaged linear polarizability (a) vs q m at 
three different frequencies. 



For small system sizes (up to about 12 sites), the 
DMRG calculations can be compared with exact results. 
This has been done for polarizability a, as well as for the 
THG coefficients (7) for various components and at dif- 
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TABLE X: Tumbling averaged THG coefficient 7 a „ in 10 3 
a.u, for different excitation frequencies and oxidation states 
of the central metal ion. 
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ferent frequencies. The accuracy is typically better than 
1% for dominant components away from resonances. 

As shown in the Fig. |S1 tumbling averaged linear polar- 
izability a av increases quadratically with increasing oxi- 
dation state, qM, on the metal ion at all the frequencies 
we have studied. However, there seems to be saturation 
in the a av value for qM > 3, at all frequencies. This 
trend seems to be consistent with the optical gaps and 
associated transition dipole moments for the low-lying 
states for different qM values, although quantitative in- 



ference about a av cannot be arrived at based on the few 
low-lying state properties we have computed. 

The tumbling averaged THG coefficients (fav), are pre- 
sented in Table [X] The THG coefficients, at all the 
three frequencies we have computed are about one or- 
der of magnitude smaller for FBP than regular metallo- 
porphines. Increase in qM leads to increase in "f av , away 
from resonances. j av value for qM = 1, 2 and 3 are 
very similar and much larger than the FBP values. For 
qM > 3, there is a large increase in j av , with qM- Thus 
metallo-porphines with large oxidation state on the metal 
ion are expected to exhibit large THG responses. 



VI. CONCLUSION 

In Summary, we have employed the symmetrized 
DMRG method for studying the linear and nonlinear 
optical properties of FBP and metallo-porphines within 
an interacting 7r-electron model. We note that the com- 
puted singlet and triplet low-lying excited state energies 
and their charge densities are in excellent agreement with 
experimental as well as many other theoretical results. 
From our bond order calculation, we conclude that por- 
phine has 18-annulenic structure in the ground state and 
the molecule gets expanded upon excitation. We have 
modelled the metalloporphine by taking an effective elec- 
tric field due to the metal ion and computed the exci- 
tation spectrum. Metalloporphines have D±h symmetry 
and hence more degenerate excited states. The ground 
state of Metalloporphines show 20-annulenic structure, 
as the charge on the metal ion increases. The linear po- 
larizability seems to increase with the charge initially and 
then saturates. The same trend is observed in third order 
polarizability coefficients. 
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